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ABSTRACT 

Context. Broad absorption line quasars (BAL QSOs) are objects showing absorption from relativistic outflows, with velocities up to 
0.2c. These manifest, in about 15% of quasars, as absorption troughs on the blue side of UV emission lines, such as Civ and Mgii. 
The launching mechanism and duration of these outflows is not clear yet. 

Aims. In this work, we complement the information collected in the cm band for our previously presented sample of radio loud BAL 
QSOs (25 objects with redshifts 1.7<z<3.6) with new observations at m and mm bands. Our aim is to verify the presence of old, 
extended radio components in the MHz range, and probe the emission of dust (linked to star formation) in the mm domain. 

Methods. We observed 5 sources from our sample, already presenting hints of low-frequency emission, with the GMRT at 235 and 
610 MHz. Other 17 sources (more than half the sample) were observed with bolometer cameras at IRAM-30m (MAMB02, 250 GHz) 
and APEX (LABOCA and SABOCA, 350 and 850 GHz, respectively). 

Results. All sources observed with the GMRT present extended emission at a scale of tens of kpc. In some cases these measurements 
allow us to identify a second component in the SED, at frequencies below 1.4 GHz, beyond the one already studied in the GHz 
domain. In the mm-band, only one source shows emission clearly ascribable to dust, detached from the synchrotron tail. Upper limits 
were obtained for the remaining targets. 

Conclusions. These findings confirm that BAL QSOs can also be present in old radio sources, or even in restarting ones, where 
favourable conditions for the outflow launching/acceleration are present. A suggestion that these outflows could be precursors of 
the jet comes from the fact that ~70% of our sample can be considered in a GigaHertz Peaked Spectrum (GPS) or Compact Steep 
Spectrum (CSS)-l-GPS phase. This would confirm the idea proposed by other authors that these outflows could be recollimated to 
form the jet. Comparing with previous works in the literature, dust emission seems to be weaker than the what expected in ‘normal’ 
QSOs (both radio loud and radio quiet ones) suggesting that a feedback mechanism could inhibit star formation in radio-loud BAL 
QSOs. 

Key words. Quasars: absorption lines - Galaxies: active - Galaxies: evolution - Radio continuum: galaxies 


1. Introduction 


In the context of Active Galactic Nuclei (AGN) outflows from 
the central region are commonly detected as absorption lines in 
more bands. In the UV/Optical range, they can be present in 
~70% of type 1 AGN, with an extension up to kpc scales and 
velocities up to ~1000 km/s ( [Harrison et al.|[2014[ ). A similar 
fraction (~60%) has been found in the quasar (QSO) population 
( [Ganguly & Brotherton||2008 l. In the X-ray band, more highly 
ionized outflows are detected as Ultra Fast Outflows (UFO) both 
in radio quiet ( Pounds et al.[2003| l and radio loud AGN ( [Tombesi] 
et al.|2014|l with much higher velocities, in the range 0.03-0.4c. 


A feedback effect from outflows on the host galaxy has been 


proven to exist by different authors in the past years (e.g. Fer- 
uglio et al.[2010[ Wang et al.[2010[ Sturm et al.[201 l[l, and lately 


they have been proven to hamper star formation (Tombesi et al.[ 
[20T5] |. 


BAL QSOs are among the objects presenting the fastest out¬ 
flows. These are detected in about 15% of QSOs as broad absorp¬ 
tion troughs in the UV spectrum, on the blue side of emission 
lines from ionised species, mainly C iv and Mg ii . They can be 


both detached or superimposed to the emission peak, and reach 
relativistic velocities of up to 0.2c ( [Hewett & Foltz|2003[ ). [All^ 
[et al.[ ( [201 l[ l have found a dependence with redshift of the BAL 
fraction, decreasing a factor of 3.5+0.4 from z~4 to ~2. The 
mechanism at the origin of these violent outflows has not been 
unveiled yet. The two main scenarios discussed in the literature 
tend to ascribe the BAL phenomenon to: 1) young objects, in 
which the strong nuclear starburst activity is still expelling a dust 
cocoon ( [Briggs et al.|1984][Sanders[2002||Farrah et al.|2007[ ), or 
2) normal QSOs, whose outflows intercept the line of sight of 
the observer ( [Elvis[[2000| . In this case, relativistic outflows are 
supposed to be commonly present in QSOs, but detected only 
when orientation is favourable. Variability of the BAL troughs 
has been explored by many authors in the past years, thanks to 
the increasing amount of available spectroscopic surveys data 


( Gibson et al.|2008] 2010[ Capellupo et al.|2()TT 2012 Vivek et 
[al.|2012[), and a typical duty cycle of about a thousand years for 


the BAL-producing outflow has been found (Filiz Ak et al. 2012 
[20T3] |. 

Several works have been published in the past years, trying 
to collect information in the different electromagnetic bands. In 
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particular, the emission in the radio band has been used to probe 
the orientation and age of these objects ([Montenegro-Montes et 


al. 2008| DiPompeo et al. 2011 Bruni et al. 2012| 20131. No 


that jets in radio loud AGN can have multiple phases of ac- 


Jamrozy 2009 Nandi et al. 2014 

with a duty cycle that de- 

pends on the source radio power 

Best et al. 

2005[ [Shabala et[ 


clear hints of a favoured scenario arose from these works, re¬ 
sulting in indications of different possible orientations and dif¬ 
ferent ages for radio loud (RL) BAL QSOs. In this work, we 
present follow-up observations of sources from our previously 
studied sample ( [Bruni et al. 2012| l. We explored the emission 
properties at m- and mm-wave lengths, to complement the multi¬ 
wavelength view of these objects, already studied at cm-band in 
our previous work. 

The detection of a strong MHz emission can be safely inter¬ 
preted as the presence of old extended radio plasma connected 
to a forme r AGN radio-active phase , that can be as old as 10^- 
10* years (Konar et al.|2006 2013 i. It has been shown indeed 


The main difference of the above samples compared to our 
target sample is the radio-loudness of our sources. With the 
radio-data presented in Bruni et al. ( 2012| l we were able to char¬ 
acterise the synchrotron spectra of our sources, and thus to study 
the peak frequency and spectral index distributions with respect 
to the ‘normal’ QSO population. Also an upper limit to the syn¬ 
chrotron emission at mm-wavelengths can be constrained. 

A search for HFP, i.e. the youngest known radio sources with 
the highest turnover frequencies, shows only a very small per¬ 


centage of sources with peak frequencies close to 20 GHz (Dal 
jlacasa et al.|2000]l with the most extreme case at 25 GHz, lead¬ 


ing to a formal age of some 50 years only (jOrienti & Dallacasaj 


2008 1 . As any upturn towards an even higher peak frequency 


would be visible in our Spectral Energy Distributions (SEDs) 
we can safely assume that the extrapolated synchrotron emis¬ 
sion reflects its true contribution at 250 GHz and that any ob¬ 
served excess emission can be attributed to the presence of cold 


ready found from literature data, for the same sample studied 
in this paper, and presented in [Bruni et al.| ( |20T^ . This kind of 
emission adds a significant information in the framework of the 
presented models for BAL QSOs. In particular, if the young- 
scenario would be the most realistic one, no further components 
than the one peaking in the MHz-GHz range should be present, 
since GigaHertz-Peaked Spectrum sources (GPS) and Compact 
Steep Spectrum sources ( CSS,[0’Dea|1998T), tog ether with High 
Erequency Peakers (HEP, Dallacasa et al.[[2()00i are among the 


youngest radio sources. To date systematic searches of diffuse 
emission around GPS and CSS sources find about 20% detec¬ 
tions (Stanghellini et al. 1990 2005[l. In light of this we use 


low-frequency GMRT observations to probe possible extended 
emission around BAL QSO with the aim of possibly confirming 
or discarding the youth scenario described in the previous para¬ 
graph. The radio phase itself does not seem to introduce signif¬ 
icant differences in RL BAL QSOs with respect to Radio Quiet 
(RQ) ones ([Bruni et al.|2014[ Rochais et al.|2014 1, thus being a 


dust. Moreover, the variability study presented in Bruni et al. 

! 2012 I excludes any possible significant variability even at high 
requencies over a 3 years time scale, for this sample of objects. 

The outline of the paper is as follows: in Sect.|^we describe 
the BAL QSO sample. The radio observations are reported in 
Sect.[^ In Sect. we present the results concerning morphology 
at MHz frequencies and dust abundance. Sect.|^is a discussion 
of the results, in the context of recent works about BAL QSOs. 

The cosmology adopted throughout the work assumes a flat 
universe and the following parameters: Hq = 71 km s ' Mpc ', 
Oa = 0.73, Qm = 0.27. 


2. The RL BAL QSO sample 

The radio-loud BAL QSO sample studied in this paper is pre¬ 
sented in [Bruni et al. ( [2012[ l. All sources were chosen among ob¬ 
jects from the 4th edition of SDSS Quasar Catalogue ( Schneid^ 
et al.|20()7 1 , drawn from the 5th data release of the SI 


tal Sky Survey (SDSS-DR5; |Adelman-McCarthy et aL 


oan Digi- 
2007] ). To 


valid tool to study general phenomenology of these objects. 

The continuum emission of dust, in the rest-frame far- 
infrared domain (EIR), can be detected at mm-wavelengths (over 
100 GHz) for objects with z~2. Objects enshrouded by gas and 
dust can host star-formation regions ( [Zahid et al.[2014[ l, and thus 
show high star-formation rates, that may indicate a young age of 
the galaxy. A flux density excess in the EIR could be an indi¬ 
cator of a different age for BAL QSOs with respect to the non- 
BAL QSO population, and thus help in discriminating between 
the orientation and the evolutionary models. There are two ma¬ 
jor works presenting (sub)mm observations on samples of BAL 
QSOs. [Willott et 3T]p003[ l showed SCUBA measurements on 
a sample of 30 radio-quiet BAL QSOs and conclude that there 
is no difference between BAL QSOs and a comparison sample 
of non-BAL QSOs. Nevertheless, [Priddey et al.[ ( [2007[ l based 
on SCUBA observations of 15 BAL QSOs found tentative ev¬ 
idence for a dependence of submm flux densities on the equiv¬ 
alent width of the characteristic C iv BAL which ‘suggests that 
the BAL phenomenon is not a simple geometric effect [... ] but 
that other variables, such as evolutionary phase, [... ] must be 
invoked’ . \CslO Orjales et al.[ ( [2012[ l discuss the far-infrared prop¬ 
erties and star-formation rates (SLR) of BAL QSOs (without dis¬ 
tinguishing between RL and RQ ones), using data from the Her- 
schel-AfTLAS project. They found no differences with respect to 
non-BAL QSOs, concluding that a scenario in which BAL QSOs 
are objects expelling a dust cocoon is improbable. 


select radio-loud objects, we cross-matched the SDSS with the 


EIRST (Eaint Images of the Radio Sky at Twenty-cm; Becker et 
[al.[2001[ l and only those with a counterpart lying <2 arc sec away 
and having S i. 4 ghz > 30 mJy were considered. All of these sat 


isfy the radio-loudness definition by Stocke et al.[ ( 19921. More¬ 
over, the selection has been limited to those objects whose red- 
shifts lie in the range 1.7 < z < 4.7 allowing the identification 
of both Civ and Mgii absorption features on SDSS spectra. In 
order to select genuine BAL QSOs, only objects with an Absorp¬ 
tion Index (Al) > 100 km/s were considered, and onW troughs 
broader than 1000 km/s were used for this calculatiorQ This re¬ 
sulted in 25 radio-loud BAL QSOs. Eor a complete description 
of the sample and the selection procedure refer to [Bruni et al.[ 

( |20T2] i. 


3. Radio observations and data reduction 

In this paper, we present observations complementary to the ones 
performed by [Bruni et al. ( 2012| i at cm wavelengths. The GMRT, 
the APEX single-dish and the IRAM 30-m telescope were used 
to extend the available SEDs extension. Table [T] summarises the 
different runs and observing setups. 


' i.e. we adopted an Al defined as Al = ^ (1 - ^)-Crfv, as in 

[et al.|200^ but with C=1 only for contiguous troughs > 1000 km/s 
null otherwise. 


Hall 
, and 
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Table 1: Summary of the observations and setups presented in this paper. IRAM and APEX telescopes made use of bolometer 
receivers. 


Run 

Date 

Telescope 

Frequency 

(GHz) 

Bandwidth 

(MHz) 

FWHM 

(arcsec) 

Number of sources 

1 

9-11 Jan 2010 

GMRT 

0.235,0.610 

33 

15-45, 4-10 

5 

2 

26 Oct-23 Nov 2010 

IRAM-30m 

250 

- 

11 

11 

3 

05-22 Aug 2010 

APEX 

850 

- 

19, 8 

4 

4 

01-09 Sep, 05 Nov 2011 

APEX 

345, 850 

- 

19, 8 

2 

5 

04-05 Jun, 13-14 Aug 2012 

APEX 

850 

- 

19, 8 

3 


3.1. Giant Metrewave Radio Teiescope 

Observations at frequencies of 235 MHz and 610 MHz with the 
GMRT were performed for five sources, during January 2010. 
We used the double frequency mode to observe simultaneously 
at the two frequencies. The total bandwidth for each band was 33 
MHz, divided into 256 channels of 0.13 MHz each. We observed 
in snapshot mode to improve the UV coverage for the sources. 
Standard phase and amplitude calibration were performed, using 
3C 286 as primary calibrator every ~4 hours and suitable phase 
calibrators near targets every ~30 minutes. Correlation was per¬ 
formed using the GSB software correlator at NCRA. Data were 
reduced with the AlPl^package, using the standard procedures. 
Flux densities were extracted from images via Gaussian fit of the 
components, using task JMFIT inside AlPS. 

3.2. IRAM-30m singie dish 

We could observe 11 sources of the BAF QSO sample at 250 
GHz with the IRAM-30m telescope, during the 2010 summer 
pool session. We used the MAMB02 117-pixel bolometer in 
ON/OFF mode, since all of our sources are point-like for this 
telescope (HPBW=11 arcsec). With average atmospheric condi¬ 
tions the detector could reach a noise of ~1 mJy/beam in ~40 
min of observing time: we observed each source for this dura¬ 
tion in order to obtain a detection or a 3-cr upper limit. Sky dip, 
calibration and pointing scans were regularly performed during 
the runs and every time the observing direction in the sky sig¬ 
nificantly changed in elevation. Focus was repeated at sunrise 
and sunset. A standard reduction was done using the MOPSlcj^ 
script provided by IRAM. 


3.3. APEX 

From 2010 to 2012 we observed with APEX a total of 9 south¬ 
ern sources from the BAF QSO sample. All of them were ob¬ 
served with the SABOCA bolometer array ( |Siringo et al.|2010) l 
at 850 GHz in photometry mode (HPBW ~8 arcsec) and two 
of the sources (0044-1-00 and 1404-1-07) also in mapping mode. 
In addition, one of the sources (0044-1-00) was observed with 
the FABOCA bolometer array ( [Siringo et al.|2009| at 345 GHz 
in photometry mode (HPBW ~19 arcsec). APEX observations 
were carried out in service mode, with typical integration times 
of 1 h per source, in order to reach RMS values around 20 
mJy/beam. Calibration was based on observations of primary 
calibrators (Mars and Uranus) as well as sky dips measured at 
the same azimuth of the targets to derive atmospheric opacity. 
A standard reduction was done using the version 2.15-1 of the 


^ http://www.aips.nrao.edu/index.shtml 
^ http://www.iram.es/IRAMES/mainWiki/CookbookMopsic 


CRUSpj^software which offers an improved pipeline for photo¬ 
metric data as compared to earlier versions. 

3.4. Error determination 

In the flux density error calculation, different contributions were 
considered for the GMRT interferometric observations: 

- The thermal noise, A5noise, which is estimated from the 
map, in empty regions of sky surrounding the target; 

- The fractional calibration error, A5caiib, estimated as the 
visibilities dispersion of the flux density calibrators; 


In particular, we followed the approach proposed by Klein et 
al. (2003). The expression used is reported below: 

AI - ,/(A5cahb ■ 5)2 + (A5noise)2 ' (1) 

V ^beam 

where Abeam and Asm are respectively the area of the synthesised 
beam and the aperture used to extract the source flux density. 
From their ratio we determine the number of beams contained in 
the source. 

For APEX and IRAM-30m data, obtained with bolometer 
receivers, the error was calculated using the respective pack¬ 
ages for data reduction, estimating the noise from off-source sub¬ 
scans. 


3.5. The poiarimetry campaign 


During 2011, we conducted a poiarimetry campaign using the 
EVFA and the Effelsberg-lOOm single dish on this same sam¬ 
ple, to implement data later presented in Bruni et al. ( 2012| l, and 
probe with deeper observations the polarisation of the faintest 
sources. The results from this campaign will be presented in a 
future paper, but we decided to use part of the obtained total flux 
density measurements to improve the SED coverage of this work 
(see Tab. Also, one of the sources observed with the EVFA 
turns out to have a resolved structure (0849H-27), and we present 
here the map (see Sect. |4.1[ ). This same source was also observed 
during our mm campaign (se e Tab.[4|). Ob servations and data re¬ 
duction were conducted as in [Bruni et al.|(|201^. 


4. Results 

In the following, we present the results from our observing cam¬ 
paign: the morphology could be studied with the GMRT inter¬ 
ferometer, while images from the IRAM and APEX bolometers 

http://www.submm.caltech.edu/ sharc/crush/ 
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1.4 GHz FIRST 



ARC SEC 

Levels = 0.5»(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam 


4.86 GHz JVLA 



ARC SEC 

Levels= 0.07»(-l, 1, 2, 4, 8, 16, 32, 64, 128, 256) mJy/beam 


5 GHz 


VLBA 



Levels = 0.2»(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam 


Fig. 1: Maps of BAL QSO 0849+27 from the FIRST survey (left panel - 1.4 GHz, beam 5.40x5.40 arcsec), our EVLA observations 
(central panel - 4.86 GHz, beam 1.95x1.25 arcsec), and our VLBA observations (right panel - 5 GHz, beam 2.99x1.15 mas). 
Contours are multiples of 3-cr, according to label. Dashed contours are negative. The synthesised beam size is shown in the lower 
left corner of each map. Note that the scale of the right panel is in mas. 


Table 2: Top lines; revised flux densities for sources presented in 
|Bruni et ar] ( |2012| l and in this work. Bottom lines: flux densities 
from our polarimetry campaign, used for this work. 


ID 

(J2000) 

Frequency 

[GHz] 

S 

[mJy] 

Telescope 

0756+37 

43 

5.2+0.8 

VLA 

0816+48 

1.4 

70.9+0.7 

VLA 

1335+02 

43 

8.4+L3 

VLA 

0842+06 

8.35 

19.5+0.8 

Effelsberg-lOOm 

0849+27 

4.86 

27.6+0.7 

EVLA 


were used for photometric measurements. The collected infor¬ 
mation is presented in Tab. while SEDs, including flux den¬ 
sities at cm wavelengths from [Bruni et al.|(|2012[ ) and our po¬ 
larimetry campaign, are presented in Fig. and discussed in 
Sect. IQ 


4.1. Morphology 

From the GMRT and EVLA maps we were able to investigate 
the morphologies of the sources at arcsec scales. The frequency 
range explored with the GMRT allowed us to put some con¬ 
straints to the presence of extended, old, radio components. 


GMRT maps 

We observed with the GMRT the five sources from our BAL 
QSO sample showing the strongest low-frequency emission in 
the flux densities collected from archival surveys data. The goal 
was the detection of extended emission at 235 or 610 MHz, in¬ 
dicating a previous radio-activity period of the central AGN, and 
thus putting a constraint on the age. Maps show components with 
deconvolved dimensions greater than zero in most cases, cor¬ 
responding to a fraction of the beam (see Eig. [^. One source 
(1159+01) presents an elongated structure at 235 MHz, con¬ 
firmed by the detection of a second component (B) at higher res¬ 
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olution in the 610 MHz map. This structure is compatible with 


the one seen in the pc-scale maps obtained by Hayashi et al. 
( |2013| l, where a jet extension at a comparable position angle is 
visible. This could confirm the presence of two different radio 
phases, at different scales, as also highlighted by the SED of this 
object (see Sec. |4.2[ ). Quantities for all sources at the two fre¬ 
quencies are presented in Tabletogether with projected linear 
sizes. 

These measurements confirm the presence of a low- 
frequency, older radio component in some BAL QSOs, thus ex¬ 
cluding that they are a subclass of young radio objects. The size 
of these components is significantly larger than the values of a 
few kpc measured for the high-frequency, unresolved, counter- 
20^, thus suggesting different emitting 


parts (see Bruni et al. 


regions for the two. The flux densities found nicely fit with col¬ 
lected data from surveys (see Eig.[^. 


EVLA map of 0849+27 


During our polarimetry campaign, we could observe the peculiar 
BAL QSO 0849+27 with the EVLA at 4.86 GHz. The map of the 
resolved structure of this source at 1.4 GHz, obtained from the 
EIRST survey, was already presented in a previous work from 
our group ( Bruni et ar||2012[ ). This resulted to be the most ex¬ 
tended BAL QSO of our sample (44 arcsec, 382 kpc, between 
components A and C). We could obtain a map at 4.8 GHz from 
our subsequent EVLA observation in 2011, taking advantage of 
the improved performance of this instrument, and also a high res¬ 
olution (pc-scale) map fo the core component from our VLBA 
program (see Bruni et ar]|2013| l. Results are presented in Eig. 

Eour out of five components detected in the EIRST data are 
visible in our EVLA map (A, B, C, E), while the flux density 
of component D seems to drop below the 3-cr significance level. 
Total flux densities as measured at 4.86 GHz are 27.6+0.7 mJy, 
2.0+0.2 mJy, 0.49+0.05 mJy, and 0.41+0.05 mJy for component 
A, B, C, E, respectively. Components E was not classified in our 
previous EIRST map, but given the clear detection we obtained 
in the EVLA map, we extracted the flux density at the corre- 
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Levels 8.S*(-1, 1, 2, 4, 8, 16, 32, 64) mjy/beam Levels 30.7*(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam Levels 15.5*(-1, 1, 2, 4, 8, 16, 32, 64) mjy/beam 



Levels 62.7*(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam Levels 5.3*(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam 



Arc seconds 

Levels 5.7*{-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam 



Arc seconds 

Levels 30.2*(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam 



60 40 20 0 -20 -40 -60 

Arc seconds 


Levels 48.3*(-l, 1, 2, 4, 8, 16, 32, 64) mjy/beam 



Arc seconds 

Levels 11.9*(-1, 1, 2, 4, 8, 16, 32, 64) mjy/beam 


Fig. 2; Maps of 5 BAL QSOs observed with the GMRT at 235 and 610 MHz. Contours are multiples of 3-cr, according to the label. 
Dashed contours are negative. The synthesised beam size is shown in the lower left corner of each map. 


spending position in the FIRST map, where a single contour was 
present. This results to be 2.36+0.15 mJy. The obtained spec¬ 
tral index for these components, are -0.54+0.04, -1.10+0.17, 
-1.78+0.18, and -1.41+0.22 for components A, B, C, and E, 
respectively. A flat spectral index (> -0.5) usually identifies the 
core for non-Doppler-boosted components: component A shows 
a spectral index compatible with that value within the error. 


while B, C, and E have a steep spectral index (< -0.5). Erom our 
VLBA observations we could confirm that core is component A, 
since at pc-scale it shows a core-jet structure, with one of the 
two components (Al) having a flat spectral index of -0.12+0.20 
between 5 and 8.4 GHz. 

The peculiar morphology of this source, explored at both kpc 
and pc scales, suggests a jet precession, or radio-activity phases 
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Table 3; Information collected with the GMRT interferometer: Col. 2 and 8 are flux densities from the maps. Col. 3-6 and 9-12 are 
deconvolved major and minor axes at 235 and 610 MHz, respectively. The projected linear size in kpc is also given. When a null 
deconvolved size were found, we considered the corresponding beam size as upper limit. For the only resolved source (1159-1-01) 
we give values for both components. 


ID 

S 

(mJy) 

Maj. axis 
(arcsec) 

235 MHz 
Min. axis 
(arcsec) 

Maj axis 
(kpc) 

Min. axis 
(kpc) 

Component 

610 MHz 

S Maj. axis 

(mJy) (arcsec) 

Min. axis 
(arcsec) 

Maj axis 
(kpc) 

Min. axis 
(kpc) 

0756+37 

<85.2 

- 

- 

- 

- 


112+4 

3.9 

2.2 

32.4 

18.3 

1159+01 

720+65 

19.6 

7.0 

164 

58.4 

A 

312+10 

4.6 

2.3 

39.5 

19.8 







B 

74.8+4.2 

10.2 

0.9 

86.3 

7.6 

1159+06 

872+81 

6.9 

<10.1 

58.7 

<85.9 


163+5 

6.7 

2.4 

57.0 

20.4 

1406+34 

249+26 

14.8 

<11.0 

122 

<90.9 


168+5 

3.9 

1.1 

32.2 

9.1 

1624+37 

104+16 

11.1 

<11.0 

84.9 

<84.1 


98.6+4.9 

4.3 

3.7 

32.9 

28.3 


with different jet axes. In fact, trajectories connecting the core 
(A) with the other components detected at kpc scale are all dif¬ 
ferent, not permitting to identify one as the counter-jet of the 
other. Moreover, the pc-scale structure shows a further launching 
direction for the most recent component (A2), not corresponding 
to any among the ones at kpc scale. In this scenario, the BAL- 
producing outflow would be present in a source presenting mul¬ 
tiple ongoing radio phases: this would not be easily explicable 
with the young scenario. 


4.2. SEDs shape from 235 MHz to 850 GHz 


We present here a study of the SED shape for the objects in 
this work. We implemented the new flux densities in the ones 
from Bruni et al. ( |2012| l, spanning from 74 MHz up to 43 GHz, 
in order to improve the overall frequency coverage. For three 
sources (0756H-37, 0816-1-48, 1335-1-02) we provide here a re¬ 
vised flux density for the VLA measurements presented in Bruni 
|et al.| ( |2OT2) l: for source 0816-1-48, during previous data reduction 
the vicinity of a strong source led to an incorrect phase refer¬ 
encing, resulting in an overestimated flux density measurement. 
While for sources 0756H-37 and 1335-1-02, our flux-extraction al¬ 
gorithm missed the source, displaced of a few arcsec from the 
map centre, due to atmospheric effects. The revised values are 
given in Tab. 


Flux densities at mm wavelengths 

Flux density from the dust grey-body thermal emission can be 
described by the following equation ([Hughes et al.|1997|): 


^obs ^ 


(1+z) 

Dl 




( 2 ) 


where is the observed flux density at a given frequency Vgi,s, 
and Vrest is the rest-frame frequency, Di is the luminosity dis¬ 
tance, z is the target redshift, is the mass absorption coeffi¬ 
cient at a given rest-frequency, B is the black-body Planck func¬ 
tion, Td is the dust temperature, and Md is the dust mass. Typical 
AGN values for Td and Md from the l atest works in the literatur e 
are 10<7’^<60 K, and Md ~ 10^Mq ( Kalfountzou et aL]|2014| l, 


while is usually calculated scaling down to the desired wave¬ 
length values estimated in previous works (e.g. k^^°^“=0.077 


Vi wavelength as a power law: 


m"^kg ‘ in Dunne et al. 2000 2011 1 , approximating the trend 


- 


4 ^ 


(3) 
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where is the dust emissivity index ( jda Cunha et al.|[2008| l. 
We could collect flux densities at mm-wavelengths for a total 
of 17 objects from the BAL QSO sample. In most cases only 
upper limits could be derived (see Tab. |^, nevertheless provid¬ 
ing useful constraints for the dust abundance. Only one source 
(0756H-37) shows a flux density at 250 GHz clearly distinguish¬ 
able from the expected synchrotron-emission tail: this is the only 
genuine example of dusty BAL QSO found among these objects. 
Despite the fact that we could not collect enough detections in 
the mm-band to perform a full fit of the dust emission using Eq. 
1^ we can still compare the obtained upper limits with the re¬ 
sults from previous works in the literature ( jOmont et al.||2003j 
[Kalfountzou et al.|2014| see Sec.[5.2[for a discussion). 


SED fitting 

For the SED fitting in the m-/cm-wavelengths domain, we 
adopted the same basic approach as in our previous work, i.e. a 
linear and a parabolic fit in logarithmic scale, determining which 
of the two functions fits best the data. The former is a simple 
model of power-law synchrotron emission, commonly used to 
fit the optically-thin part of the emission above the peak, for old 
sources peaking at frequencies below the available ones. The lat¬ 
ter is a first-order approximation of an optically-thick (on the 
left side of the peak) plus an optically-thin synchrotron emis¬ 
sion (on the right side of the peak) capable of fitting the SED 
of a young radio source, peaking in the MHz-GHz range. Upper 
limits were excluded from the dataset used for fitting, as well 
as datapoints clearly belonging to a second component in the 
MHz range (if any). After introducing the new data, we obtained 
best fits with the same functions (line or parabola) than previ¬ 
ously adopted, except for sources 0842-1-06, that resulted to have 
a parabolic shape, and 1335-1-02, that after the introduction of the 
revised flux density from Tab. show a HEP component peak¬ 
ing at ~20 GHz (see Fig. [^. Another source (1159-1-06) showed 
a second synchrotron component in the MHz range, peaking at 
about 100 MHz, while 0816H-48, that was not fitted in [Bruni ^ 
[ar] ( [2012| l, now presents a parabolic component in the GHz range 
using again the revised flux density at 1.4 GHz. 

In two of the three sources detected at 250 GHz the emis¬ 
sion is most probably due to synchrotron component (1237H-47, 
1406-1-34), since it perfectly fits in the shape of the emission tail. 
For source 0756-1-37 we found a significant excess of emission at 
250 GHz (2.0+0.5 mJy), with respect to the expected contribu¬ 
tion of the synchrotron emission at the same frequency. Thus we 
can consider this emission as most probably produced by dust. 
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Table 4: Results for the 17 BAL QSOs observed in the mm-band: 
250 GHz flux densities from IRAM-30m, 345 and 850 GHz from 
APEX. We give 3-cr upper limits in case of non-detection. 


Name 

S250 

(mJy) 

S 345 
(mJy) 

S 850 

(mJy) 

0044+00 

- 

<260 

<30 

0756+37 

2.0+0.5 

- 

- 

0816+48 

<1.8 

- 

- 

0842+06 

- 

- 

<71 

0849+27 

<2.4 

- 

- 

1014+05 

<3.0 

- 

- 

1102+11 

<3.0 

- 

- 

1159+01 

- 

- 

<150 

1159+06 

<6.0 

- 

<161 

1229+09 

<3.3 

- 

<66 

1237+47 

4.6+1.0 

- 

- 

1304+13 

<3.0 

- 

- 

1327+03 

<3.6 

- 

<201 

1335+02 

- 

- 

<78 

1337-02 

- 

- 

<66 

1404+07 

- 

- 

<90 

1406+34 

9.3+0.8 

- 

- 


5. Discussion 

In the following, we discuss our findings and put them in the 
context of the works present in the literature. 


5.1. Low-frequency components 

A consistent number of restarting radio sources have been found 
in the past 20 years, and hypotheses about their nature have been 
proposed ( [Czerny et al.|200^ Wu|2009 Marecki et al.|2006 1 . A 
fraction corresponds to GPS sources or CSS sources, known to 
be young radio sources. In some cases, up to three radio phases 
are detectable in the same object ( Brocksopp et al.||2007 1. |Cz^ 
jerny et ak] ( |2009| ) associate the intermittent activity of the cen¬ 
tral engine with the radiation pressure instability of the accre¬ 
tion disk. Concerning our results in the MHz range, the fact that 
hints of old, extended, radio components were found, despite the 
limited number of objects observed with the GMRT, could in¬ 
dicate that an age of 10^-10® years for the BAL-hosting QSOs 
is not rare. This finding is in line with what arose from our pre¬ 
vious works ( [Bruni et al.||2012| [2013[ l, where we discussed the 
presence of old components in the SED of our sample. In some 
cases a radio-restarting scenario, and the complex dynamics in¬ 
volved by that, could be invoked to explain a two-component 
SED (1159 -hOI, 1159-h06, 1335-1-02, 1406h- 34), showing both a 
GHz-peaked and a MHz component. In other cases (1014-1-04, 
1229-1-09, 1304-1-13, 1327-1-03) a peak in the MHz range could 
still be present, but not seen because of insufficient frequency 
coverage. Considering this, and the SEDs of objects previously 
studied in [Bruni et al. ( 2012[ l, a ~70% of our sample could be in 
a GPS or GPS-hCSS phase. 

Once again this shows how BAL-producing outflows can 
be present not only in young radio sources, but also in more 
complex scenarios. Eor example, sources found to have multi¬ 
ple radio phases ongoing (CSS-hGPS), should have gone through 
an unstable radiation pressure phase, causing an intermittent 
BAL-producing outflow acceleration (favoured during the high- 
pressure phases). This could include young, just-started radio 
sources (GPS/HEP) as well as restarted radio sources (e.g. those 


with a CSS-hGPS SED). If this link with ignition/recollimation 
of the radio jet would be confirmed, an outflow collimation to 
form a jet (as already proposed by Elvis|2000) l could be invoked 
to explain the BAL variability. 


5.2. Dust emission 

Several works in the past years, especially with the advent of 
the Herschel space observatory, has constrained the dust emis¬ 
sion properties in AGN, and verified the connection with star- 
formation rate. Here we provide a comparison between our re¬ 
sults for BAL QSOs in the mm-band and what found by other 
authors. 

A very similar study to ours, in terms of instrumentation, sen¬ 
sitivity, and setup, is the one presented by Omont et al. ( 2003[ l. 
They performed 250-GHz observations of 35 optically luminous 
RQ QSOs {Mb < -27.0), with a similar redshift range to ours 
(1.8< z <2.8), performed with the MAMBO bolometer at the 
same frequency. They found that 26+9% of the sources present 
an emission at that frequency. Since they reached an RMS very 
similar to our observations (~1 mJy) we can compare our de¬ 
tection rate with this percentage: only 1 out of 17 (~6%) of our 
sources shows 250-GHz emission attributable to dust, a substan¬ 
tially smaller fraction than the one found by these authors. 

More recently, Kalfountzou et al. ( 2014| published results of 
the Herschel-ALTLAS project regarding EIR properties of radio- 
loud and radio-quiet QSOs. Using five different bands (100, 160, 
250, 350, 500 pm) they fitted the dust emission for both groups. 
At the mean redshift of our sample (z~2.3) they found a mean 
flux density for radio-loud QSOs of 23.5+2.1 mJy at 350 pm 
(corresponding to our 850 GHz observations), and 21.3+2.8 mJy 
at 500 pm (~600 GHz). These are values below the sensitivity we 
reached with APEX, but using a simple power law to extract the 
expected mean flux density at 250 GHz from the previous values, 
we obtain a value of ~ 17 mJy, that would be well detectable with 
the RMS of ~1 mJy we could reach at IRAM-30m. This could 
suggest that our sample of RE BAL QSOs is also poorer in dust 
content than normal RE QSOs, in addition to RQ ones (as seen 
comparing with Omont et al.|2003| l. 

In the light of these works, our results support the idea of 
BAL QSOs not being specially dusty objects. In addition, al¬ 
though with modest statistical significance, our work suggests 
that dust emission in the RE BAL QSO sub-class could be 
weaker than expected. [Tombesi et al.[ ( [2()15[ l found that fast out¬ 
flows (~0.25c) from the accretion disk in AGN can hamper star 
formation, impacting the interstellar medium. BAL-producing 
outflows, although at a lower ionization stage than the ones con¬ 
sidered by those authors, can present velocities up to 0.2c, and, 
in the light of these results, could have a similar effect on the 
SLR of the host galaxy. This needs to be investigated with fur¬ 
ther observations, that could confirm this on larger samples of 
RE BAL QSOs. 


6. Conclusions 

We performed observations in the m-band with the GMRT of 
5 RE BAL QSOs, that already showed hints of emission in the 
MHz range, and in the mm-bands for 17 RE BAL QSOs, from 
our previously studied sample. We aimed at exploring the emis¬ 
sion in the low-frequency regime, and the grey-body emission 
from dust, respectively. The conclusions are the following: 

- All 5 objects observed at low frequencies present emission 
from extended components, indicating the presence of an 
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old radio emission. In some cases a restarting radio-activity 
can be invoked to explain the double component, MHz and 
GHz-peaked, present in their SEDs. This could suggest an 
intermittent BAL phase, associated with periods of radio¬ 
restarting activity. This is also supported by the morphology 
of some objects, both from this work (0849-1-27) and from 
|Bruni et al.| ( |2013| l, and by the fact that ~7Q% of our RL 
BAL QSOs sample can be considered in a GPS or CSS-hGPS 
phase, thanks to the data presented here. 

- Only 1 out of 17 sources (~6%) presents a clear contribu¬ 
tion at 250 GHz from the dust grey-body emission. In the 
other cases 3-cr upper limits have been derived. Comparing 
our results with the fraction of dust-rich RQ QSOs found by 
Omont et al.| ( |2003| ), resulting in a percentage of ~26%, we 
found that RL BAL QSOs do not present a larger fraction 
of dust-rich objects with respect to the RQ QSO population. 
Also comparing with more recent works from the Herschel- 
ATLAS collaboration ( jKalfountzou et aL]|2014[ ), we found 
a lack of dust emission with respect to mean values for RL 
QSOs. Since the amount of dust can be connected to star- 
formation and thus to the age of the host galaxy, this results 
suggest that RL BAL QSOs are not a hosted by particularly 
young galaxies. The fact that, despite of radio loudness, these 
objects present even less dust emission (and consequently a 
lower star formation rate) than RQ QSOs, could suggest that 
BAL-producing outflows are able to hamper star formation 
in the host galaxy. 

- Both the obtained results, from observations performed at 
m- and mm-wavelengths, suggest that BAL QSOs are not 
commonly young radio objects, or objects still expelling 
their dust cocoon from the central region. They could 
be radio-restarting objects, which present relativistic out¬ 
flows in conjunction with some periods of favourable emis¬ 
sion/acceleration conditions. 
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Frequency (GHz) 


Fig. 3: SEDs of the 18 BAL QSOs observed during the m/mm-wavelengths campaign (x-axis: GHz; y-axis: mJy). 235 MHz and 
610 MHz flux densities from the GMRT; 250 GHz flux densities from IRAM-30m; 345 and 850 GHz flux densities from APEX. 
Measurements at other frequencies are taken from |Bruni et al.| ( 2012[ ), or from our polarimetry campaign (see Sect. 3). Triangles are 
3-cr upper limits. Solid lines are parabolic or linear fits, according to the criteria discussed in Sec. 4.2. 
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